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QO (5ZLAb*tf*acl; An electromagnetic survey method for surveying an area of sea floor that is thought or known to contain a subter- 
O ranean hydrocarbon reservoir comprising obtaining a first survey data set with a vertical electric dipole (VED) antenna for generating 
^ vertical current loops and a second survey data set with a vertical magnetic dipole (VMD) antenna for generating horizontal current 
^ loops. In an alternative embodiment, the VMD antenna is dispensed with and the horizontal electromagnetic field is derived from the 
naturally occurring magnctolelluric (MT) electromagnetic field. In another alternative embodiment, the VED data is compared with 
a background geological model instead of VMD or MT data. The invention also relates to a survey apparatus comprising VED and 
VMD antennae, to planning a survey using this method, and to analysis of survey data taken using this survey method. The first and 
second survey data sets allow the galvanic contribution to the detector signals collected at a detector to be independently contrasted 
£^ with the inductive effects. This is important to the success of using electromagnetic surveying for identifying hydrocarbon reserves 
^* and distinguishing them from other classes of structure. 
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TITLE OF THE INVENTION 

ELECTROMAGNETICSUKVEYIN rT FOR HYD ROCARBON RESERVOIRS ^^ 
BACKGROUND OF THE INVENTION 

The invention relates to seafloor electromagnetic surveying for oil and other 
5 hydrocarbon reserves. 

Determining the response of the sub-surface strata within the earth's crust to 
electromagnetic fields is a valuable tool in the field of geophysical research. The 
geological structures associated with thermally, hydrothermally, technically or 
magmatically active regions can be studied. In addition, electromagnetic surveying, or 
10 sounding, techniques can provide valuable insights into the nature, and particularly the 
likely hydrocarbon content, of subterranean reservoirs in the context of subterranean 
oil exploration and surveying. 

Seismic techniques are often used during oil exploration to identify the 
existence, location and extent of reservoirs in subterranean rock strata. Whilst seismic 
15 surveying is able to identify such structures, the technique is often unable to 
distinguish between the different possible compositions of pore fluids within them, 
especially for pore fluids which have similar mechanical properties. In the field of oil 
exploration, it is necessary to determine whether a previously identified reservoir 
contains oil or just aqueous pore fluids. To do this, an exploratory well is drilled to 
20 determine the contents of the reservoir. However, this is an expensive process, and 
one which provides no guarantee of reward. 

Whilst oil-filled and water-filled reservoirs are mechanically similar, they do 
possess significantly different electrical properties and these provide for the possibility 
of electromagnetic based discrimination testing. A known technique for 
25 electromagnetic probing of subterranean rock strata is the passive magneto-telluric 
(MT) method. The signal measured by a surface-based electromagnetic detector in 
response to electromagnetic (EM) fields generated naturally, such as within the earth's 
upper atmosphere, can provide details about the surrounding subterranean rock strata. 



WO 2004/008183 




PCT/GB2003/002379 



-2- 



However, for deep-sea surveys, all but those MT signals with periods corresponding 
to several cycles per hour are screened from the seafloor by the highly conductive 
seawater! Whilst the long wavelength signals which do penetrate to the seafloor can 
be used for large scale undersea probing, they do not provide sufficient spatial 
5 resolution to examine the electrical properties of the typically relatively small scale 
subterranean reservoirs. Moreover, since MT surveying relies primarily on 
horizontally polarised EM fields, it is intrinsically insensitive to thin resistive layers. 

Nonetheless, measurements of electrical resistivity beneath the seafloor have 
traditionally played a crucial role in hydrocarbon exploration and reservoir assessment 

10 and development. In industry, subterranean resistivity data have generally been 
obtained almost exclusively by wire-line logging of wells. There are, though, clear 
advantages to developing non-invasive geophysical methods capable of providing 
such information from the surface or seafloor. Although inevitably such methods 
would be unable to provide comparable vertical resolution to wireline logging, the 

15 vast saving in terms of avoiding the costs of drilling test wells into structures that do 
not contain economically recoverable amounts of hydrocarbon would represent a 
major economic advantage. 

In research fields that are not of commercial interest, geophysical methods for 
mapping subterranean resistivity variations by various forms of EM surveying have 

20 been in use for many years [1, 2, 3, 10]. Proposals for finding hydrocarbon reservoirs 
using such EM surveying have also been made [4, 5] and applications to the direct 
detection of hydrocarbons using horizontal electric dipole (HED) sources and 
detectors have proved successful [6, 7]. 

To successfully map subterranean resistivity variations in the field of oil 

25 exploration, the orientation of the current flows induced by EM signals must be 
carefully considered [6]. The response of seawater and subterranean strata (which will 
typically comprise planar horizontal layers) to EM signals is generally very different 
for horizontally and vertically flowing current components. For horizontally flowing 
current components, the coupling between the layers comprising the subterranean 
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strata is largely inductive. This means the presence of thin resistive layers (which are 
indicative of hydrocarbon reservoirs) do not significantly affect the EM fields detected 
at the surface since the large scale current flow pattern is not affected by the thin layer. 
On the other hand, for vertical current flow components, the coupling between layers 
5 is largely galvanic (i.e. due to the direct transfer of charge). In these cases even a thin 
resistive layer strongly affects the EM fields detected at the surface since the large 
scale current flow pattern is interrupted by the resistive layer. It is known therefore 
that significant vertical components of induced current are required to satisfactorily 
perform an EM survey in the field of oil exploration. 

10 However, sole reliance on the sensitivity of vertical current flow components 

to the presence of a thin resistive layer can lead to ambiguities. The effects on detected 
EM fields arising from the presence a thin resistive layer can be indistinguishable 
from the effects arising from other realistic large scale subterranean strata 
configurations. In order to resolve these ambiguities, it is known that it is necessary to 

15 determine the response of the subterranean strata to both horizontal (i.e. inductively 
coupled) and vertical (i.e. vertically coupled) induced current flows [6]. 

Hence it is important when designing a practical EM survey for detecting 
buried hydrocarbon reservoirs to distinguish between source and detector 
configurations in which the coupling between layers is largely inductive due to 

20 horizontal currents (in which case the survey has little sensitivity to the presence of a 
thin reservoir) and those in which the coupling between layers is largely galvanic due 
to vertical currents (in which case blocking of the passage of this current flow by a 
reservoir leads to a survey which is strongly sensitive to the presence and boundary of 
hydrocarbon within the reservoir). 

25 Figure la schematically shows a surface vessel 14 undertaking EM surveying 

of a subterranean strata configuration according to a previously proposed method [6]. 
The subterranean strata configuration includes an overburden layer 8, an underbidden 
layer 9 and a hydrocarbon layer (or reservoir) 12. The surface vessel 14 floats on the 
surface 2 of the seawater 4. A deep-towed submersible vehicle 19 carrying a HED 
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antenna 21 is attached to the surface vessel 14 by an umbilical cable 16 providing an 
electrical and mechanical connection between the deep-towed submersible vehicle 19 
and the surface vessel 14. The HED antenna broadcasts a HED EM signal into the 
seawater 4. 

5 One or more remote detectors 25 are located on the seafloor 6. Each detector 

25 includes an instrument packages 26, a detector antenna 24, a floatation device 28 
and a ballast weight (not shown). In practice, each detector antenna 24 will generally 
comprise an array of antenna elements, for example, a pair of orthogonal dipole 
antennae elements. The detector antenna 24 measures a signal in response to EM 

10 fields induced by the HED antenna in the vicinity of the detector 25. The instrument 
package 26 records the signals for later analysis. 

The HED antenna 21 generates both inductive and galvanic current flow 
modes with the relative strength of each mode depending on HED antenna-detector 
geometry. At detector locations which are broadside to the HED antenna axis, the 

1 5 inductive mode dominates the response. At detector locations which are in-line with 
the HED antenna axis, the galvanic mode is stronger [6, 8, 9, 10]. The response at 
detector locations in both the in-line and broadside configurations is governed by a 
combination of the inductively and galvanically coupled modes and these tend to work 
in opposition. 

20 Figure lb shows in plan view an example survey geometry according to the 

previously proposed method in which sixteen detectors 25 are laid out in a square grid 
on a section of seafloor 6 above a subterranean reservoir 56 having a boundary 
indicated by a heavy line 58. The orientation of the subterranean reservoir is indicated 
by the cardinal compass points (marked N, E, S and W for North, East, South and 

25 West respectively) indicated in the upper right of the figure. To perform a survey, a 
source starts from location c A' and is towed along a path indicated by the broken line 
60 through location until it reaches location 'C which marks the end of the survey 
path. As is evident, the tow path first covers four parallel paths aligned with the 
North-South direction to drive over the four "columns" of the detectors. This part of 
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the survey path moves from location 'A' to 'B\ Starting from location C B\ the survey 
path then covers four paths aligned with the East- West direction which drive over the 
four "rows" of detectors. Each detector is thus driven over in two orthogonal 
directions. The survey is completed when the source reaches the location marked *C\ 
5 During the towing process, each of the detectors 25 presents severed different 

orientation geometries with respect to the source. For example, when the source is 
directly above the detector position Dl and on the North-South aligned section of the 
tow path, the detectors at positions D5, D6 and D7 are at different ranges in an end-on 
position, the detectors at positions D2, D3 and D4 are at different ranges in a 
10 broadside position and the detector at positions D8 and D9 are midway between. 
However, when the source later passes over the detector position Dl when on the 
East- West aligned section of the tow path, the detectors at positions D5, D6 and D7 
are now in a broadside position, and the detectors at position D2, D3 and D4 are in an 
end-on position. Thus, in the course of a survey, and in conjunction with the 

15 positional information of the source, data from the detectors can be used to provide 
details of the signal transmission through the subterranean strata for a comprehensive 
range of distances and orientations between source and detector, each with varying 
galvanic and inductive contributions to the signal propagation. In this way a simple 
continuous towing of the source can provide a detailed survey which covers the extent 

20 of the subterranean reservoir. 

This previously proposed method has been demonstrated to provide good 
results in practice. However, some limitations of the method have been identified. 

Firstly, since the two modes cannot be easily separated there will generally be 
a level of cross-talk between them at a detector and this can lead to ambiguities in the 

25 results. 

Secondly, in order to obtain survey data from both in-line and broadside 
geometries, the HED antenna needs to be re-oriented at each HED antenna survey 
location. This requires the surface vessel to make multiple passes over broadcast 
locations and can lead to complex and long tow patterns. 
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Thirdly, a HED antenna based EM survey can only provide the best data 
possible at discrete survey locations. This is because of the geometric requirements of 
a HED antenna survey which dictate that, at any point during the survey, data can only 
be optimally collected from those detectors to which the HED antenna is arranged 
5 either in-line or broadside. At other orientations, separation of the inductively and 
galvanically coupled signals becomes more much difficult and data are less reliable. 
For instance, referring to Figure lb, when the HED antenna is at a point on the tow 
path above the detector marked Dl and on the North-South aligned section of the tow 
path, in-line data can only be collected from the detectors marked D5, D6 and D7, 

10 whilst broadside data can only be collected form the detectors marked D2, D3 and D4. 
The other detectors provide only marginally useful information at this point of the 
survey. Furthermore, if the HED antenna is at the location identified by reference 
numeral 57 in Figure lb, which is on a North-South aligned section of the tow path, 
in-line data can be collected from the detectors marked D3, D8, D9 and D10, but 

1 5 broadside data cannot be collected from any of the detectors. Since both broadside and 
in-line data are required for optimal analysis, the best data possible with the square 
detector array shown in Figure lb can only be collected from points along the tow 
path directly above the detector locations. 
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SUMMARY OF THE INVENTION 

According to a first aspect of the invention there is provided an 
electromagnetic survey method for surveying an area that potentially contains a 
5 subterranean hydrocarbon reservoir, comprising: obtaining a first survey data set using 
a vertical electromagnetic source signal; and obtaining a second survey data set using 
a horizontal electromagnetic source signal. 

In this case, references to vertical and horizontal indicate that a significant 
component, preferably a major component, of the respective signals should be aligned 

10 with the vertical and horizontal axes. It is not necessary that the signals are perfectly, 
or even closely, aligned to the vertical and horizontal axes, although fairly close 
alignment is preferred to provide a strong signal and reduce the complexity of 
analysis, for example alignment within +/-30 0 is desirable. 

The vertical electromagnetic source signal is preferably provided by a vertical 

15 electric dipole (VED) source. The VED source can be towed over the survey area 
above the seafloor by a surface vessel and interconnecting umbilical cable. 

The horizontal electromagnetic source signal is preferably provided by a 
vertical magnetic dipole (VMD) source. The VMD source can also be towed over the 
survey area above the seafloor by a surface vessel and interconnecting umbilical cable. 

20 Most preferably, the VED source and VMD source are part of a single 

submersible tow vehicle and the collection of the first and second survey data sets is 
done concurrently, typically by towing the vehicle over the seafloor along a predefined 
survey path. During surveying, the tow vehicle can be switched at short intervals 
between activation of the VED source and the VMD source. Alternatively, the VED 

25 and VMD sources can be operated simultaneously at different frequencies to allow 
easy separation of the two signals at the detectors. In either case, the complete survey 
can be carried out with a single tow, i.e. both the first and second survey data sets can 
be obtained concurrently. In contrast, in previously proposed EM survey methods 
based on use of a survey vehicle with a HED antenna, two sets of survey data had to 
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be obtained with different tows over different paths in order to probe the response of 
an area of interest galvanically and inductively. Being able to collect the first and f 
second survey data sets concurrently is of great practical advantage, since it obviates 
the need to perform two separate tows over different survey paths. 
5 As an alternative to using a VMD source, the horizontal electromagnetic 

source signal could be provided by magnetotelluric (MT) electromagnetic fields, 
which are naturally occurring fields in the earth's upper atmosphere. In this case, the 
tow vehicle need only include a VED source. (A VED source for determining the 
resistivity and extent of large scale sedimentary structure is known [3, 16]. The known 
10 VED source extends over the full depth of the seawater in which it operates and is 
driven at the quasi-static limit.) If the second survey data set is obtained from MT 
electromagnetic fields, the data can be obtained concurrently with collection of the 
first survey data set or, more likely, at a different time. 

Since a VED antenna generates electromagnetic field configurations which are 
15 sensitive to the 'presence of a thin resistive layer, whereas a VMD generates 
electromagnetic field configurations which are sensitive to the background structure, 
by operating the VED and VMD antennae at different frequencies, or consecutively at 
the same frequency, the galvanically and inductively coupled modes of current 
transfer can be easily separated, so avoiding the difficulty of the 'mixed-mode' fields 
20 generated in previously proposed survey methods. 

Furthermore, the vertical electromagnetic signal induces significant 
components of electric current which are normal to a surface of a body of water in 
which the electromagnetic survey is performed This limits the deleterious effects of 
the finite seawater depth compared to those seen in previously proposed 
25 electromagnetic survey methods. The vertical electromagnetic signal can therefore be 
operated in shallower water than has previously been possible and allows 
electromagnetic surveying to be performed in areas which have not previously been 
amenable to electromagnetic surveying. 
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If the first survey data set is obtained by towing a source of the vertical 
electromagnetic signal along a tow-path relative to an array of detectors within the 
area of interest, a comprehensive survey can be performed rapidly. Furthermore, 
because patterns of current flow induced by vertical electric dipole (and vertical 
5 magnetic dipole) electromagnetic signals are cylindrically symmetric, detector data 
can be consistently and reliably collected from all detectors within the array for all 
source locations. This has not been possible with previously proposed survey methods 
for which data can only be reliably collected for specific relative orientations of the 
source and detector. Furthermore, if the layers comprising the subterranean strata are 
10 largely homogeneous and isotropic, the circular symmetry of the vertical 
electromagnetic signal can obviate the need to determine the orientation of the 
detectors within the area of interest. In cases where the layers comprising the 
subterranean strata are not homogeneous and isotropic, knowledge of the detector 
antenna orientation, for instance measured by a compass, will allow more thorough 
1 5 data analysis to be performed. 

According to a second aspect of the invention, there is provided a method of 
analysing results from an electromagnetic survey of an area potentially containing a 
subterranean hydrocarbon reservoir, comprising: providing a first survey data set 
obtained using a vertical electromagnetic source signal; providing a second survey 
20 data set obtained using a horizontal electromagnetic source signal; generating a first 
normalisation data set specific to the first survey data set; generating a second 
normalisation data set specific to the second survey data set; combining the first 
survey data set and first normalisation data set to obtain a first results data set that 
represents a difference between the first survey data set and the first normalisation 
25 data set; and combining the second-survey data set and second-normalisation data set 
to obtain a second results data set that represents a difference between the second 
survey data set and the second normalisation data set. 

The survey method distinguishes between the presence of structure within the 
subterranean strata which influence galvanically coupled electric currents flowing 
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vertically between layers (for example thin resistive layers which can be indicative of 
the presence of hydrocarbon), and the presence of structures which influence 
inductively coupled currents flowing horizontally within layers. Differences between 
the first survey data set and the first normalisation data set reflect structures in the 
5 subterranean strata which influence both galvanically and inductively coupled electric 
currents, whereas differences between the second survey data set and the second 
normalisation data set reflect the presence only of structures which influence 
inductively coupled currents. Comparisons between the first and second data sets 
therefore allow the two classes of structure to be distinguished and hence the presence 

10 of thin resistive layers (such as hydrocarbons) to be inferred. Interpretation of the first 
and second survey data sets may be improved if they are normalised to the respective 
first and second normalisation data sets or first and second functions specific to the 
first and second data sets respectively. 

The vertical electromagnetic source signal is preferably a vertical electric 

1 5 dipole (VED) electromagnetic signal and the horizontal electromagnetic source signal 
is preferably a vertical magnetic dipole (VMD) electromagnetic signal. Alternatively, 
the horizontal electromagnetic source signal could be a magnetotelluric (MT) 
electromagnetic field signal. 

The first results data set may represent the difference between the first survey 

20 data set and the first normalisation data set as a function of position within the area 
and the analysis of the first results data set may include identifying a location of a 
boundary of the subterranean hydrocarbon reservoir 

The normalisation data sets or functions may be calculated from a rock 
formation model or from the first and second survey data sets. 

25 According to a third aspect of the invention there is provided a computer 

program product bearing machine readable instructions for implementing a method of 
analysing results from an electromagnetic survey as described above. 



WO 2004/008183 




PCT/GB2003/002379 



-11- 



According to a fourth aspect of the invention there is provided a computer 
apparatus loaded with machine readable instructions for implementing the method of 
analysing results from an electromagnetic survey as described above. 

According to a fifth aspect of the invention there is provided a method of 
5 planning an electromagnetic survey of an area that potentially contains a subterranean 
hydrocarbon reservoir, comprising: creating a model of the area to be surveyed, 
including a seafloor, a rock formation containing a hydrocarbon reservoir beneath the 
seafloor and a body of water above the seafloor; setting values for water depth, depth 
below the seafloor of the hydrocarbon reservoir, and resistivity structure of the rock 
10 formation; performing a simulation of an electromagnetic survey in the model of the 
survey area by calculating first and second survey data sets for a vertical 
electromagnetic source signal and a horizontal electromagnetic source signal. 

The vertical electromagnetic source signal is preferably a vertical electric 
dipole (VED) electromagnetic signal and the horizontal electromagnetic source signal 
15 is preferably a vertical magnetic dipole (VMD) electromagnetic signal. Alternatively, 
the horizontal electromagnetic source signal could be a magnetotelluric (MT) 
electromagnetic field signal. 

Repeated simulations for a number of source-to-detector distances and 
frequencies can be performed in order to allow optimum surveying conditions in 
20 terms of source-to-detector distance and frequency of EM signal for probing the 
hydrocarbon reservoir to be selected when performing an electromagnetic survey. The 
effects of differing detectors array configurations and source tow paths can also be 
modelled. The model for the simulation may include a body of air above the body of 
water, wherein the simulation takes account of signal propagation paths including the 
25 body of air when calculating the first and second survey data sets. The first and second 
survey data sets may be normalised relative to respective first and second 
normalisation data sets or functions specific to the first and second survey data sets 
respectively. 
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According to a sixth aspect of the invention there is provided a computer 
program product bearing machine readable instructions for implementing the method 
of planning an electromagnetic survey as described above. 

According to a seventh aspect of the invention there is provided a computer 
5 apparatus loaded with machine readable instructions for implementing the method of 
planning an electromagnetic survey as described above. 

According to an eighth aspect of the invention there is provided a submersible 
vehicle for subsea electromagnetic surveying, comprising: a vertical electric dipole 
antenna; and a vertical magnetic dipole antenna, 
10 The vertical and magnetic electric dipole antennae are preferably mounted 

such that respective dipole axes thereof are aligned. 

In the case that the antennae are driven by signals generated by a power supply 
in the surface vessel, the submersible vehicle will typically include a waveform 
generator such as a cycloconverter operable to transform a high voltage, low current 
15 AC drive signal received from an umbilical cable into a low voltage, high current AC 
drive signal to drive the VMD and VED antennae. 

According to a ninth aspect of the invention there is provided a survey 
apparatus comprising: a submersible vehicle according to the eighth aspect of the 
invention; a signal power supply unit for generating a high voltage, low current AC 
20 signal power supply for the submersible vehicle; and an umbilical cable releasably 
connectable at ends thereof to the signal power supply unit and the submersible 
vehicle. The survey apparatus can be assembled ready for use by connecting the 
umbilical cable to the signal power supply unit, which will typically be located on the 
surface vessel, and the submersible vehicle, which can then be deployed into the sea. 
25 According to a tenth aspect of the invention there is provided a surface vessel 

carrying a survey apparatus according to the ninth aspect of the invention. 

According to an eleventh aspect of the invention there is provided an 
electromagnetic survey method for surveying an area that potentially contains a 
subterranean hydrocarbon reservoir, comprising: obtaining a first survey data set using 



/ 
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a first electromagnetic field that drives electric current in vertical planes; and 
obtaining a second survey data set using a second electromagnetic field that drives 
electric current in horizontal planes. 

According to a twelfth aspect of the invention there is provided a method of 
5 analysing results from an electromagnetic survey of an area that is thought or known 
to contain a subterranean hydrocarbon reservoir, comprising: providing a survey data 
set obtained from a vertical electric dipole (VED) source; generating a normalisation 
data set specific to the survey data set, wherein the normalisation data set is calculated 
from a rock formation model; and combining the survey data set and normalisation 

10 data set to obtain a results data set that represents a difference between the survey data 
set and the normalisation data set. In this aspect of the invention, the normalisation 
data set is of sufficient quality as regards taking account of spatial variation in 
resistivity (in 1, 2 or 3 dimensions) to obviate the need to compare the VED data with 
VMD or MT data. The rock formation model preferably includes resistivity, and may 

1 5 be derived from a combination of geological data and resistivity data. The geological 
data can be from seismological surveying and the resistivity data from well logging. 
Other sources, such as neutron data could also be used. 

According to a thirteenth aspect of the invention there is provided a method of 
planning an electromagnetic survey of an area that is thought or known to contain a 

20 subterranean hydrocarbon reservoir, comprising: creating a model of the area to be 
surveyed including a seafloor, a rock formation containing a postulated hydrocarbon 
reservoir beneath the seafloor, and a body of water above the seafloor; setting values 
for water depth, depth below the seafloor of the postulated hydrocarbon reservoir, and 
resistivity structure of the rock formation; and performing a simulation of an 

25 electromagnetic survey in the model of the survey area by obtaining a survey data set 
from a simulated vertical electric (VED) dipole source; and adjusting the model to 
remove the postulated hydrocarbon reservoir and repeating the simulation to obtain a 
normalisation data set for comparison with the survey data set. 
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According to a fourteenth aspect of the invention there is provided a method of 
monitoring an area that contains a subterranean hydrocarbon reservoir by 
electromagnetic surveying, comprising: obtaining a first VED survey data set of the 
area; waiting a period of time; obtaining a second VED survey data set; and 
5 combining the first and second VED survey data sets to obtain a results data set that 
represents a difference between the first and second VED survey data sets, said 
difference being indicative of changes of the reservoir during said period of time. 

Optionally, the first and second VED survey data sets can also combined with 
a normalisation data set. Preferably, each of the first and second VED survey data sets 

10 are independently normalised before combining them, typically using the same 
normalisation data set, although different normalisation sets could be used. It would 
also be possible to normalise only one of the two data sets. Another possibility would 
be to obtain the results data set before normalisation and then normalise the combined 
results data set. In other cases, no normalisation may be necessary. The normalisation 

15 data set can be obtained from a survey of the area by a vertical magnetic dipole 
(VMD) source or obtained from a survey of the area by naturally occurring 
magnetotelluric (MT) electromagnetic fields or calculated from a rock formation 
model, for example. 

Typically the period of time between VED surveys will be at least a few days 

20 and usually of the order of weeks to a few months. For example, the period of time 
between obtaining VED survey data sets for comparison can be at least a day, a week, 
or a month. 

The method is particularly useful for identifying changes in boundaries of the 
reservoir during the period of time through analysis of the changes in the VED survey 
25 data, for example to monitor reservoir depletion. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of the invention and to show how the same may be 
carried into effect reference is now made by way of example to the accompanying 
5 drawings. 

Figure la shows in schematic vertical section a surface vessel undertaking an 
EM survey according to a previously proposed method; 

Figure lb is a schematic plan view showing an example survey geometry 
according to the previously proposed method in which sixteen detectors are laid out 
10 on a section of seafloor above a subterranean reservoir; 

Figure 2a shows in schematic vertical section a surface vessel undertaking an 
EM survey according to an embodiment of the invention; 

Figure 2b shows a vertical electric dipole (VED) antenna and a vertical 
magnetic dipole (VMD) antenna for use in the EM survey shown in Figure 2a; 
15 Figure 3a schematically shows some instantaneous characteristic electric and 

magnetic field lines associated with a VED antenna; 

Figure 3 b schematically shows some instantaneous characteristic electric and 
magnetic field lines associated with a VMD antenna; 

Figure 4 schematically shows in vertical section a model subterranean strata 
20 configuration; 

Figure 5a is a graph schematically showing the strength of the vertical and 
radial components of electric field seen at a detector in response to a VED EM signal 
for two different subterranean strata configurations; 

Figure 5b is a graph schematically showing the phase of the vertical and radial 
25 components of electric field seen at a detector in response to a VED EM signal for 
two different subterranean strata configurations; 

Figure 6a is a graph schematically showing the strength of the azimuthal 
component of electric field seen at a detector in response to a VMD EM signal for two 
different subterranean strata configurations; 
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Figure 6b is a graph schematically showing the phase of the azimuthal 
component of electric field seen at a detector in response to a VMD EM signal for two 
different subterranean strata configurations; 

Figure 7 is a graph schematically showing the normalised field strengths of the 
5 vertical, radial and azimuthal components of the electric fields shown in Figures 5a 
and 6a; 

Figure 8 schematically shows in vertical section another model subterranean 
strata configuration; 

Figure 9a is a graph schematically showing the strength of the vertical and 
10 radial components of electric field seen at a detector in response to a VED EM signal 
for two different subterranean strata configurations; 

Figure 9b is a graph schematically showing the phase of the vertical and radial 
components of electric field seen at a detector in response to a VED EM signal for 
two different subterranean strata configurations; 
15 Figure 10a is a graph schematically showing the strength of the azimuthal 

component of electric field seen at a detector in response to a VMD EM signal for two 
different subterranean strata configurations; 

Figure 10b is a graph schematically showing the phase of the azimuthal 
component of electric field seen at a detector in response to a VMD EM signal for two 
20 different subterranean strata configurations; 

Figure 11 is a schematic plan view showing an arrangement of sixteen 
detectors on a section of seafloor above a subterranean reservoir, an example EM 
source tow path is also shown. 

Figure 12a is a graph showing the normalised strength of the vertical 
25 component of electric field seen at a detector in response to a VMD EM signal for 
different depths of seawater; 

Figure 12b is a graph showing the normalised strength of the radial component 
of electric field seen at a detector in response to a VMD EM signal for different 
depths of seawater; 
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Figure 13a schematically shows in vertical section another model subterranean 
strata configuration; 

Figure 13b schematically shows in horizontal section the model subterranean 
strata configuration shown in Figure 13a; and 
5 Figure 14 is a graph schematically showing the strength of the vertical and 

radial components of electric field seen at a detector in response to a VED EM signal 
for three different subterranean strata configurations. 
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PF t An irn INSCRIPTION 

Figure 2a of the accompanying drawings schematically shows a surface vessel 
14 undertaking EM surveying of a subterranean strata configuration (or rock 
5 formation) according to an embodiment of the invention. Features shown in Figure 2a 
which are functionally similar to features shown in Figure la are given the same 
reference numeral but shall be described again for the sake of completeness. The 
subterranean strata configuration includes an overburden layer 8, an underburden layer 
9 and a hydrocarbon layer (or reservoir) 12. The surface vessel 14 floats on the surface 
10 2 of the seawater 4. A deep-towed submersible vehicle 19 carrying an EM source 18 is 
attached to the surface vessel 14 by an umbilical cable 16 providing a detachable 
electrical and mechanical connection between the deep-towed submersible vehicle 19 
and the surface vessel 14. The deep-towed vehicle 19 includes an echo location 
package 27 which assists in maintaining the EM source 18 at an appropriate height 
15 above the seafloor 6. The umbilical cable may be a standard deepwater remotely 
operated vehicle cable including three high-voltage conductors and three optic fibres. 
The EM source 18 includes a vertical electric dipole (VED) antenna 22 for 
broadcasting a first electromagnetic source signal in the form of a VED EM source 
signal (giving rise to VED EM signals within the seawater and subterranean strata 
20 configuration) and a vertical magnetic dipole (VMD) antenna 23 for broadcasting a 
second electromagnetic source signal in the form of a VMD EM source signal (giving 
rise to VMD EM signals within the seawater and subterranean strata configuration). 
The EM source 18 receives electrical power from the ship's on-board power supply 

via the umbilical cable 16. 
25 A signal power supply unit 30 on the surface vessel 14 provides signal power 

supplies to a waveform generator in the form of a cycloconverter 29 via the umbilical 
cable 16. In this example, the cycloconverter 29 is in the deep-towed submersible 
vehicle. The cycloconverter 29 generates suitable VED and VMD alternating current 
(AC) drive signals from the signal power supplies and these are supplied to the VED 
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antenna 22 and VMD antenna 23 within the EM source 18. The VED and VMD AC 
drive signals respectively comprise a first and a second set of frequency components. 
The AC drive signals drive the VED antenna 22 and VMD antenna 23 to broadcast 
the VED and VMD EM signals into the seawater 4. In this example, the VED antenna 
5 and the VMD antenna are driven (via the cycloconverter) by separate signal power 
supplies from the signal power supply unit 30. The separate signal power supplies are 
carried to the cycloconverter 29 on one of the high voltage conductors within the 
umbilical cable and are returned together on the remaining conductor acting as a 
common earth. The signal power supply unit 30 is powered from a standard 3 -phase 

10 power supply on the surface vessel. In other examples, a signal power supply unit may 
be located within the deep-towed submersible and be supplied via the umbilical cable 
with a standard 3 -phase power supply on the surface vessel. The signal power supply 
unit generates stable-frequency high-voltage (order 2kV) AC signal power supplies 
which are transformed and switched at the deep-towed submersible vehicle by the 

15 cycloconverter. The cycloconverter includes a transformer and switching bridge to 
produce suitable high current VED and VMD AC drive signals to apply to the VED 
and VMD antennae. By driving the VED and VMD antennae at different frequencies, 
the VED and VMD EM detector signals can be easily separated. Alternatively only 
one of the VED and VMD antennae could be driven at any one time, with a routing 

20 switch being used to switch between the two antennae, for example at regular 
intervals of a few minutes. 

Another design option is to generate independent VED and VMD signal power 
supplies from separate signal power supply units housed in the surface vessel, the 
signal power supplies being supplied separately down a common umbilical cable. The 

25 VED and VMD AC drive signals could then be handled by separate cycloconverters 
(or other waveform generators) in the submersible vehicle 19, thereby providing two 
independent electrical systems for the two signals. 

* One or more remote detectors 25 are located on the seafloor 6. Each detector 
25 includes an instrument package 26, a detector antenna 24, a floatation device 28 
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and a ballast weight (not shown). The detector antenna 24 produces a VED detector 
signal and a VMD detector signal in response to the VED and VMD EM signals in the 
vicinity of the detector 25 (the terms VED and VMD detector signal are used to 
indicate which of the VED and VMD EM signals the detector signal corresponds to, 
5 and not to reflect the EM field pattern in the vicinity of the detector 25). The 
instrument package 26 records the VED and VMD detector signals for later analysis. 
The detector signals may include multiple components of the EM field (both magnetic 
and electric) in the vicinity of the detector. The detector antenna 24 in this example 
includes three orthogonal electric dipole detector antennae, a first one arranged 
10 vertically for measuring the vertical electric field component of the VED and VMD 
EM signals at the detector, and second and third orthogonal ones arranged in a 
horizontal plane for measuring orthogonal horizontal electric field component of the 
VED and VMD EM signals at the detector. The detector antenna 24 detects (and the 
instrument package 26 accordingly records) both the amplitude and phase of the 
15 electric field arisirig from the VED and VMD EM signals at the detector 25 resolved 
along these three orthogonal directions. 

Figure 2b is a schematic perspective view further detailing the EM source 18 
and the deep-towed submersible vehicle 19. The deep-towed submersible vehicle is 
attached to the umbilical cable 16 by a releasable mechanical coupling 35. As noted 
20 above, the EM source includes a VED antenna 22 and a VMD antenna 23. The VED 
antenna is attached to the deep-towed submersible vehicle by a releasable mechanical 
coupling 31. The VMD antenna is mounted on a support frame 36 which is attached to 
a lower end of the VED antenna. 

The VED antenna includes an upper VED antenna electrode 32 and a lower 
25 VED antenna electrode 34 separated by a VED antenna main body 33. The VED AC 
drive signal from the cycloconverter is coupled to the VED antenna electrodes 32, 34 
via cabling within the deep-towed submersible vehicle and VED antenna main body 
33. The conducting seawater 4 presents an unscreened return path for the electrical 
current comprising the VED AC drive signal and so provides the VED EM source 
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signal. In this example, the VED antenna 22 has a length (i.e. distance between the 
upper and lower antenna electrodes) of 100 m and is supported such that the lower 
electrode is maintained by the echo location package 27 at a height of 10 m above the 
seafloor. 

5 The VMD antenna 23 is a loop antenna and is coupled to the VMD AC drive 

signal from the cycloconverter via cabling within the deep-towed submersible vehicle, 
the VED antenna main body 33 and a length of VMD antenna connection cable 37. 
The VMD antenna may comprise a single loop, but in this example, several conductor 
turns are employed to provide a larger magnetic dipole moment for a given supply 

10 current and loop diameter. 

Whilst not shown in this example, a faired cowling may be placed around the 
VED and/or VMD antenna to reduce drag and help maintain the vertical alignment of 
the upper and lower electrodes. Stabiliser fins, again not shown in this example, may 
also be employed to assist in maintaining the stability and orientation of the VED 

15 and/or VMD antennae during towing. A short streamer or drogue may also be attached 
to the VMD to assist towing stability. At low towing speeds (for example <~0.8 ms" 1 
(<~1.5 kts) over the ground) the umbilical 16 between the surface vessel 14 and the 
mechanical coupling 35 shown in Figure 2b is approximately vertical close to the EM 
source 18. Nonetheless, standard ultra-short base line acoustic navigation and location 

20 methods can be used to determine the position of upper and lower ends of the VED 
antenna to better than a few meters and any deviation of the VED antenna from 
vertical can therefore be accounted for in subsequent data processing. It is therefore 
not essential that the VED antenna induces only vertical current flow loops, although 
the VED detector signal will be strongest, and the signal data more easily 

25 interpretable, if a major component of the current flow induced by the VED antenna is 
in a vertical plane. 

Figure 3 a is a schematic diagram showing some instantaneous characteristic 
electric and magnetic field lines associated with the VED antenna when driven by the 
VED AC drive signal from the cycloconverter. The electric field is schematically 
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shown by solid lines and marked E^ 0 , and the magnetic field is schematically shown 
by broken lines and marked B™. The electric E^ 0 and magnetic fields are 
both cylindrically symmetric about the long-axis of the VED antenna 22. The electric 
field E^ 0 drives current (i.e. has components) only along radial and vertical 
5 directions in a cylindrical co-ordinate system centred on and aligned with the dipole 
antenna 22, there is no azimuthal component of the electric field E^ 0 . The magnetic 
field B* 80 associated with the VED antenna, on the other hand, has only azimuthal 
components. Whilst not shown for simplicity, the magnetic field is not confined 
to the horizontal plane containing the centre of the VED antenna 22. An 
1 0 electromagnetic field configuration of the form shown in Figure 3a is often referred to 
as a toroidal magnetic (TM) field configuration. 

Figure 3b is a schematic diagram showing some instantaneous characteristic 
electric and magnetic field lines associated with the VMD antenna 23 when driven by 
the VMD AC drive signal from the cycloconverter. The electric field is schematically 
1 5 shown by solid lines and marked E VMD , and the magnetic field is schematically shown 
by broken lines and marked B™ 0 . The electric E™ 15 and magnetic B VMD fields are 
both cylindrically symmetric about a vector centred on and perpendicular to the plane 
of the VMD antenna loop. The electric field E VMD drives current (i.e. has components) 
only along azimuthal directions in a cylindrical co-ordinate system centred on and 
20 perpendicular to the VED antenna loop, there are no radial or vertical components of 
the electric field E™ 0 , accordingly the VMD EM source signal induces current flow 
in horizontal planes. The magnetic field B^ associated with the VMD antenna, on 
the other hand, has both radial and vertical components, but no azimuthal 
components. Whilst not shown for simplicity, the electric field E VMD is not confined 
25 to the horizontal plane containing the VMD antenna 23. An electromagnetic field 
configuration of the form shown in Figure 3b is often referred to as a toroidal electric 
(TE) field configuration. 

Since a VED generates only TM mode fields, which are sensitive to the 
presence of a thin resistive layer, whereas a VMD generates only TE mode fields, 



/ 



WO 2004/008183 




PCT/GB2003/002379 



which are sensitive to the background structure, by operating the sources at two 
slightly different frequencies, or consecutively at the same frequency, the two modes 
can be easily separated, so avoiding the difficulty of the * mixed-mode' fields 
generated in previously proposed survey methods. 
5 Whilst the electric field configuration for the VED antenna shown in Figure 3a 

includes horizontal radial components in addition to vertical components, it is the 
vertical components of electric field (which are not present in the electric field 
configuration for the VMD antenna shown in Figure 3b) which play a crucial role in 
the galvanic coupling mechanism. Accordingly, the electric field configuration for the 

10 VED antenna is sometimes referred to as a vertical electric field, although it will be 
understood that the electric field comprises vertical loops which include non-zero 
radial horizontal components 

During an EM survey, the VED and VMD EM source signals generate VED 
and VMD EM signals that propagate outwards both into the seawater 4 and 

15 downwards through the seafloor 6 and into the subterranean strata 8, 9, 12 comprising 
the subterranean strata configuration. In each case, at practical frequencies for this 
method and given the typical resistivity of the media 4, 8, 9, 12, propagation occurs by 
diffusion of EM fields. The rate of decay in amplitude and the phase shift of the VED 
and VMD EM signals are controlled both by geometric spreading and by skin depth 

20 effects. Because in general the subterranean strata 8, 9, 12 are more resistive than the 
seawater 4, skin depths in the subterranean strata 8, 9, 12 are longer. As a result, the 
detector signals seen at the seafloor 6 by the detector 25 (at EM source-detector 
separations if of greater than a few hundred metres) are dominated by the components 
of the VED and VMD EM signals which have propagated downwards through the 

25 seafloor 6, along within the subterranean strata 8, 9, 12, and back up to the detector 
25. 

At the end of the EM surveying experiment, a remotely operable release 
system allows the instrument package 26 to be detached from a ballast weight (not 
shown) so that an in-built flotation device 28 can carry the instrument package 26 to 
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the surface 2 for recovery and retrieval of data for analysis and interpretation. Both the 
amplitude and the phase of the detector signals recorded during a survey depend on 
the resistivity structure of the subterranean strata 8, 9, 12 - and so, in principal, a 
survey consisting of many source (transmitter) and detector (receiver) locations can 
5 provide a multi-dimensional image, by geophysical inversion, of subterranean 
resistivity. 

The technique described here exploits the large resistivity contrast that exists 
between a hydrocarbon reservoir (typically tens of Qm or higher) and the over- and 
under-lying sediments (typically ~ 2 Qm or less). Such a contrast has a detectable 

10 influence on controlled source electromagnetic (CSEM) data collected at the seafloor 
6 above the hydrocarbon reservoir 2. The effect of the reservoir is most detectable in 
CSEM data at an appropriate frequency, and if the horizontal range from the EM 
source 18 to the detector antenna 24 is of the order of 2 to 5 times the depth of burial 
of the reservoir 12 in typical situations. 

15 The following text describes specific EM source configurations, geometries 

and data reduction approaches that allow the effect of a hydrocarbon reservoir on the 
outcome of a controlled source electromagnetic survey to be detected and analysed in 
practice. 

20 A specific approach to combining contrasting djpole sources when p ro specting for 
hydrocarbon filled reservoirs 

The survey method shown in Figure 2a exploits the diverse properties of 
electromagnetic induction outlined above to provide a first survey data set (derived 
25 from the VED detector signals) which is sensitive to the presence of thin resistive 
layers (exploiting the largely galvanic effects associated with vertical components of 
the induced current flow) while simultaneously obtaining a second survey data set 
(derived from the VMD detector signals) which is more sensitive to the background 
sedimentary structure, but relatively insensitive to a thin hydrocarbon layer or 
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reservoir (exploiting the dominantly inductive effects associated with horizontal 
components of the induced current flow). By collecting both data types potential 
ambiguities in interpretation can be avoided. 

Unlike a single HED source (such as shown in Figure la and discussed above) 
5 which excites both the inductive and galvanic current transfer modes simultaneously, 
the EM source 18 described above excites each mode independently. The VED 22 
generates vertical current loops, and excites the galvanically coupled mode whilst the 
VMD 23 generates horizontal current loops and excites the inductively coupled mode. 
Since the VED 22 and VMD 23 can be driven independently (either at different times 

10 from each other or simultaneously) the effects of the two current transfer modes (i.e. 
galvanic and inductive) can be easily separated during data analysis. 

Figure 4 shows in schematic vertical section a first, or hydrocarbon, model 
subterranean strata configuration which is used to examine the differing responses of 
the subterranean strata to the VED and VMD EM signals. A VED antenna 22, a VMD 

15 antenna 23 and a detector antenna 24 are shown, other components of the survey 
equipment previously seen in Figure 2a are not included for simplicity, but will be 
understood from the description above. In the hydrocarbon model subterranean strata 
configuration shown in Figure 4, a section of seafloor 6 lies beneath an infinite depth 
of seawater 4 which has resistivity 0.3 Qm. The subterranean strata beneath the 

20 seafloor 6 comprise a 1 km thick overburden layer 8, a 100 m thick hydrocarbon layer 
12, and an infinitely thick underburden layer 9. The overburden layer 8 and the 
underbidden layer 9 have relatively low resistivities, in this example 1 Qm, due to 
aqueous saturation of pore spaces. The hydrocarbon reservoir 12 has a relatively high 
resistivity, in this example 100 Qm, due to the presence of the poorly-conducting 

25 hydrocarbon occupying much of the pore spaces. 

Figure 5a is a graph schematically showing the logarithm of the modelled 
radial and vertical electric field component amplitudes, Logio(E), seen by a detector in 
response to a VED EM source signal as a function of separation, or range, R, between 
the VED antenna and the detector. In this example, the VED AC drive signal is a 
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quasi-square wave at a frequency of 0.25 Hz and the electric fields are calculated per 
unit source electric dipole moment. The VED antenna is 100 m long and supported 
such that the lower electrode is 10 m above the seafloor. The model field strengths are 
shown both for the hydrocarbon model strata configuration shown in Figure 4 (solid 
5 lines) and for a second, or normalising, model subterranean strata configuration 
(dashed lines). The normalising model subterranean strata configuration is similar to 
the hydrocarbon model subterranean strata configuration seen in Figure 4, but without 
the hydrocarbon layer (i.e. a continuous infinite half-space of resistivity 1 Qm beneath 
an infinite depth of seawater with resistivity 0.3 Qm). In each case, the curves relating 
10 to the vertical and radial components of the detected field are marked as ver and rad. 
respectively. 

Figure 5b is a graph schematically showing the phase, <|>, relative to the VED 
AC drive signal used to drive the VED antenna, of the modelled radial and vertical 
electric field components plotted in Figure 5a. Again, curves are shown both for the 
1 5 hydrocarbon (solid lines) and the normalising (dashed lines) model subterranean strata 
configurations and marked ver and rad for the vertical and radial electric field 
components in the VEM detector signal as appropriate. 

It can be seen from Figure 5a that at short ranges (<~2000 m) the modelled 
radial electric fields for the hydrocarbon model subterranean strata configuration show 
20 a reduction in magnitude relative to the response of the uniform half-space 
normalising model subterranean strata configuration. This is caused by the 
interruption of current flow by the resistive hydrocarbon layer. At longer ranges 
however, the modelled radial electric field strength within the hydrocarbon model 
subterranean strata configuration is enhanced relative to the response of the 
25 normalising model subterranean strata configuration. This effect is characteristic in 
the galvanically coupled mode of a structure in which the resistivity first increases and 
then decreases with depth below the seafloor. The vertical electric field components at 
the seafloor show a similar pattern, although the decrease in the magnitude of the 
fields at shorter ranges is more pronounced and extends to around 6000 m. This is 
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because of the dependence of the vertical fields on the vertical component of current 
flow. 

It can be seen from Figure 5b that the phase of the detected signals is also 
strongly dependent on the presence of the hydrocarbon reservoir. For example, beyond 
5 about 4000 m, the rate of change of phase with increasing separation is significantly 
lower for both the radial and vertical electric field components in the hydrocarbon 
model subterranean strata configuration (solid lines) compared to the normalising 
model subterranean strata configuration (dashed lines). 

Figure 6a is a graph schematically showing the logarithm of the modelled 

10 azimuthal electric field component amplitude, Logio(E), seen by a detector in response 
to a VMD EM source signal as a function of separation, or range, R, between the 
VMD antenna and the detector. In this example, the VMD AC drive signal is a quasi- 
square wave at a frequency of 0.25 Hz and the electric fields are calculated per unit 
source magnetic dipole moment. The VMD antenna is supported 10 m above the 

15 seafloor. The modelled field strength is again shown for both the hydrocarbon (solid 
lines) and the normalising (dashed lines) model subterranean strata configurations 
discussed above. 

Figure 6b is a graph schematically showing the phase, relative to the VMD 
AC drive signal used to drive the VMD antenna, of the modelled azimuthal electric 
20 field components plotted in Figure 6a. Again, curves are shown for both the 
hydrocarbon (solid lines) and the normalising (dashed lines) model subterranean strata 
configurations. 

It can be seen from Figures 6a and 6b that the modelled field strength and the 
modelled phase are very similar for both the hydrocarbon and the normalising model 
25 subterranean strata configurations. The hydrocarbon layer has little effect oil either of 
the magnitude or phase of the azimuthal detector signal seen in response to the VMD 
EM signal (which induces horizontal current flows). 

As can be seen in Figures 5a and 6a, the modelled electric field component 
amplitudes drop rapidly with increasing source-detector separation R. In a typical EM 



WO 2004/008183 




PCT/GB2003/002379 



-28- 

survey the detector signal amplitude is therefore likely to vary by several orders of 
magnitude over a useful range of source-detector separations. The phase <|> of the 
detector signal increases steadily with increasing separation. In order to visualise more 
clearly the effects of a subterranean strata configuration on the detector signal during a 
5 survey, it is convenient to normalise the VED and VMD detector signals to first and 
second normalisation data sets respectively. The first and second normalisation data 
sets may reflect calculated detector signals for model subterranean strata 
configurations or functions specific to the VED and VMD EM source signals 
respectively. For example, one simple reference model to use would be the 
10 normalising model subterranean strata configuration described above. In other 
examples it may be preferred to use different models, for example, one similar to the 
normalising model subterranean strata configuration but including the true (rather than 
infinite) depth of seawater at the survey site. If appropriate, and if the relevant a priori 
information is available, a more complex reference model may be used, although it is 

15 generally desirable to use the simplest reference model that can broadly represent the 
large scale background properties of the subterranean strata. 

Figure 7 is a graph schematically showing the same electric field strengths 
plotted in Figures 5a and 6a for the hydrocarbon model subterranean strata 
configuration but after normalisation by the normalising model subterranean strata 

20 configuration. The normalised field strengths N are marked VED vcr and VED rad for the 
vertical and radial electric fields seen in response to a VED EM source signal (solid 
lines in Figure 5a) and VMDaz for the azimuthal field seen in response to a VMD EM 
source signal (solid line in Figure 6a). The normalisation is carried out by dividing the 
electric fields comprising the detector signals calculated for the hydrocarbon model 

25 subterranean strata configuration by those calculated for the normalising model 
subterranean strata configuration with the same source-detector geometry. Whilst not 
shown, in the case of the phase data, normalisation is carried out by subtracting the 
phase calculated for the normalising model subterranean strata configuration model 
from that of the hydrocarbon model subterranean strata configuration. 
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As an alternative, normalisation could be based on the survey data itself, for 
example using data collected adjacent to the target or applying a smoothing algorithm, 
such as a box-car mean or a low-frequency spatial filter, to the first and second data 
sets. 

5 As with Figure 5a, it can be seen from Figure 7 that the presence of the 

hydrocarbon layer leads to increased radial and vertical electric field strengths seen at 
a detector in response to a VED EM signal beyond around 4000 m (radial) and 6000 
m (vertical). For example, at around 8000 m the radial and vertical components are 
enhanced by factors of around 15 and 5 respectively. As with Figure 6a, it can also be 

10 seen from Figure 7 that the presence of the hydrocarbon layer has little effect on the 
azinnithal electric field strengths seen by a detector in response to a VMD EM signal. 
This is apparent since the normalised electric field amplitude is close to unity for all 
source-detector separations. 

Even though the VMD detector signals are not strongly affected by the 

15 presence of the hydrocarbon layer, observing the response of a subterranean strata to 
horizontally induced current flow is an important part of the method since, as noted 
above and explained further below, characterising the background structure in this 
way helps to resolve possible ambiguities in the interpretation of VED detector signal 
data alone. This is especially useful in areas of exploration where the background 

20 structure is not well understood. 

Figure 8 shows in schematic vertical section a third, or increasing-resistivity, 
model subterranean strata configuration used to highlight some of the ambiguities 
which may arise if only VED detector signal data are collected. A VED antenna 22, a 
VMD antenna 23 and a detector antenna 24 are shown. As with the hydrocarbon 

25 model subterranean strata configuration, the increasing-resistivity model subterranean 
strata configuration shown in Figure 9 includes a section of seafloor 6 beneath an 
infinite depth of seawater 4 which has resistivity 0.3 Qm. In the increasing-resistivity 
model subterranean strata configuration, the strata beneath the seafloor 6 comprise a 
500 m thick first layer 40 with resistivity 1 Qm, a 1km thick second layer 42 with 
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resistivity 5 Qm and an infinitely thick third layer 44 with resistivity 10 Qxn. The 
increasing-resistivity model subterranean strata configuration approximates to a 
subterranean strata configuration displaying a steadily increasing resistivity with 
depth. This is not an unreasonable feature of submarine sedimentary basins due to the 
5 progressive expulsion of conductive pore fluids with increasing depths by a rising 
overburden pressure. 

Figure 9a is a graph schematically showing the logarithm of the modelled 
radial and vertical electric field component amplitudes, Logio(E), seen by a detector in 
response to a VED EM source signal as a function of separation, or range, R 9 between 

10 the VED antenna and the detector. In this example, the VED AC drive signal is a 
quasi-square wave at a frequency of 0.25 Hz and the electric fields are calculated per 
unit source electric dipole moment. The VED antenna is 100 m long and supported 
such that the lower electrode is 10 m above the seafloor. The model field strengths are 
shown both for the hydrocarbon model strata configuration shown in Figure 4 (solid 

15 lines) and for the increasing-resistivity model subterranean strata configuration shown 
in Figure 8 (dashed lines). In each case, the curves relating to the vertical and radial 
components of the detected field are again marked as ver and rad respectively. 

Figure 9b is a graph schematically showing the phase, <J>, relative to the VED 
AC drive signal used to drive the VED antenna, of the modelled radial and vertical 

20 electric field components plotted in Figure 9a. Again, curves are shown both for the 
hydrocarbon (solid lines) and the increasing-resistivity (dashed lines) model 
subterranean strata configurations and marked ver and rad for the vertical and radial 
electric field components of the VED detector signal as appropriate. 

It can be seen from Figures 9a and 9b that the VED detector signals (both 

25 radial and vertical) are similar in amplitude and phase for both the hydrocarbon and 
the increasing-resistivity model subterranean strata configurations. Accordingly, 
whilst the presence of a hydrocarbon layer within an otherwise uniform subterranean 
strata configuration can be identified by observing the response of the subterranean 
strata to a VED EM source signal (as shown by the curves in Figures 5a and 5b), it 



WO 2004/008183 




PCT/GB2003/002379 



-31- 



can be difficult to distinguish between a hydrocarbon layer within an otherwise 
uniform subterranean strata configuration and a non-hydrocarbon-containing 
subterranean strata configuration displaying an increasing resistivity with depth such 
as is associated with some submarine sedimentary basins. 
5 Figure 10a is a graph schematically showing the logarithm of the modelled 

azimuthal electric field component amplitude, Logio(E), seen by a detector in response 
to a VMD EM source signal as a function of separation, or range, R, between the 
VMD antenna and the detector. In this example, the VMD AC drive signal is a quasi- 
square wave at a frequency of 0.25 Hz and the electric fields are calculated per unit 
10 source magnetic dipole moment. The VMD antenna is supported 10 m above the 
seafloor. As with Figure 9a, the modelled field strength is shown for both the 
hydrocarbon (solid line) and the increasing-resistivity (dashed line) model 
subterranean strata configurations. 

Figure 10b is a graph schematically showing the phase, <j>, relative to the VMD 
15 AC drive signal used to drive the VMD antenna, of the modelled azimuthal electric 
field components plotted in Figure 10a. Again, curves are shown for both the 
hydrocarbon (solid line) and the increasing-resistivity (dashed line) model 
subterranean strata configurations. 

It can be seen from Figure 10a that the detected azimuthal electric field 
20 amplitude seen in response to a VMD EM source signal for the hydrocarbon and the 
increasing-resistivity model subterranean strata configurations beyond a source- 
detector separation R of about 4000 m begin to differ. This allows the ambiguity 
between the two models in response to VED EM source signals alone to be resolved. 

In summary, the VED detector signals can be used to distinguish between a 
25 structured subterranean strata configuration (e.g. as shown in Figures 4 and 7) and a 
non-structured subterranean strata configuration (e.g. the uniform model subterranean 
strata configuration described above) and the VMD detector signals can be used to 
distinguish between different structured subterranean strata configurations to identify 
the presence of a hydrocarbon layer or reservoir. 
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Whilst in the above example a VMD EM signal is employed to determine the 
response of the subterranean strata configuration to horizontally induced currents, 
similar information on background structure could be gathered using marine 
magnetotelluric (MT) data [12]. As noted above, the ionospheric EM fields used as 
5 the source in MT sounding excite predominantly horizontal current flow in the earth. 
MT data are therefore notoriously insensitive to thin resistive layers, however they can 
be used to determine the background structure [7]. It is likely that in many 
circumstances VMD detector signal data of the type discussed above will be preferred 
since marine MT data are rarely collected at frequencies high enough to resolve 

10 structure at the same scale allowed when using controlled source VMD EM source 
signals. Furthermore, the effect of distant coastlines can often be seen in MT data and 
this can increase the complexity of the interpretation even if the local structure is 
simple. A further advantage of employing controlled source VMD EM source signals 
in addition to controlled source VED EM source signals is that provided the VMD 

15 antenna is driven kt a similar frequency as the VED antenna, the VMD EM signals 
provide information about background structure at the same scale as that resolved 
with the VED EM signals. 

Use of a mobile source and multiple fixed detectors 

20 

In order to perform a thorough survey over a large area, and to provide curves 
similar to those shown in Figures 5a, 5b, 6a, 6b, 8a, 8b, 9a and 9b, controlled source 
electromagnetic measurements will typically be made with many EM source and 
detector locations. Whilst these could be made in parallel using a plurality of EM 
25 sources operating at different frequencies and a plurality of detectors, or in series 
using a single EM source and a single detector which are repositioned between 
measurements, it will generally be more efficient to employ a plurality of detectors in 
fixed positions and a single repositionable EM source, as indicated in Figure 2a. 
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The EM source 1 8 shown in Figure 2a can require significant power to drive it, 
of the order tens of kilowatts, or greater for signals detectable at ranges of several 
kilometres. The umbilical cable 16 connecting the EM source 18 to the survey vessel 
14 supplies this power and makes it relatively straightforward to make the EM source 
5 mobile- It can then be towed in an appropriate survey pattern by the surface survey 
vessel 14. Since in many situations surveys of this kind are liable to take place over 
areas of the seafloor 6 where sensitive engineering installations exist or are planned, 
there are significant advantages to using a source which does not physically come into 
contact with the seafloor 6. Provided that the separation between the EM source 18 

10 and the seafloor 6 is small compared to a skin depth of the investigating field in 
seawater, the survey can still be completed. (Although it is noted that coupling of the 
VED EM signal into the subterranean strata is weaker than that of a HED [13]). 

Figure 1 1 is a schematic plan view showing an example layout of sixteen 
detectors 25 distributed across a section of seafloor 6 above subterranean reservoir 56. 

15 The reservoir 56 has a linear extent on the order of several km and its boundary is 
indicated by a heavy line 58. The orientation of the subterranean reservoir is indicated 
by the cardinal compass points (marked N, E, S and W for North, East, South and 
West respectively) indicated in the upper right of the figure. In this example, the 
detectors 25 are uniformly distributed in a square-grid pattern so as to approximately 

20 cover the reservoir 56. In performing a survey, an EM source (not shown) starts from 
the location marked 'A' in Figure 1 1 and is towed, whilst broadcasting continuously 
as described above, along a path indicated by the broken line 60, the survey is 
completed when the source reaches the location marked *B\ VED and VMD detector 
signal data are continuously collected by the detectors 25 throughout the towing 

25 process and the position of the EM source relative to the network of detectors is also 
logged. 

During the towing process, each of the detectors 25 presents several different 
source-detector separations and orientations relative to the EM source. Accordingly, 
by following the tow path marked, EM survey detector signal data are collected for 
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many different source-detector separations along many different directions. These 
detector signal data can be inverted to provide a thorough map of the subterranean 
strata configuration using appropriate geophysical inversion techniques. In this way a 
simple continuous towing of the source can provide a detailed survey which covers 
5 the extent of the subterranean reservoir 56. 

Figure lb shows for comparison an example tow path followed by a previously 
proposed survey geometry using a HED antenna in an electromagnetic survey of the 
same area as shown in Figure 11. Features of Figure lb similar to those shown in 
Figure 1 1 will be understood from the above. However, in this example, the tow path 

10 is longer. Whilst again stating at the location marked 'A', the tow path shown in 
Figure lb is approximately twice as long, ending at the location marked 4 C\ as that 
shown in Figure 11. When performing an EM survey according to the present 
invention, detector signal data can be collected from all detectors when the EM source 
is at all points along the tow path shown in Figure 1 1, irrespective of orientation. 

15 However, the geometric requirements of a HED antenna survey are such that 

at any point during the survey, data can only be reliably collected from those detectors 
to which the HED antenna is arranged either in-line or broadside, as described further 
above. Accordingly, with a HED antenna survey, not only is a longer tow path 
required for a given area of interest but much less data can be usefully collected along 

20 this path. It is also clear that broadside and in-line data which are to be analysed 
together (i.e. relating to the same source and detector locations) cannot be collected 
simultaneously. 

Although the above example is based on a square detector grid, it will be 
understood that a wide variety of detector placements may be used. For example other 
25 high symmetry regular grids, such as triangular or rectangular, may be used. In 
addition irregular grids may be used that have no high level of symmetry. 

In the case of the detectors, there are further advantages in using static devices. 
Firstly, since detectors suitable for the EM surveying are generally internally powered 
and relatively cheap compared to a suitable EM source, a plurality of detectors can 
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easily by arrayed around an area of interest so as to provide multiple source-detector 
ranges and positions for a single EM source position as shown in Figure 11. Secondly, 
the task of the detectors is to measure electric or magnetic components of the VED 
and VMD EM signals at the seafloor 6. In typical applications, the signal-to-noise 
5 ratio of the VED and VMD detector signals is important to the success and resolution 
of the survey, and should be maximised. Moving a detector inevitably generates noise, 
whether the field components comprising the detector signals are magnetic or electric. 
In the case of electric fields, any motion of the detector through the conducting 
seawater 4 in the presence of earth's geomagnetic field generates an electromotive 
10 force across the detector antenna 24. Detector movements will therefore map into 
spurious detector signals in the recorded data. In the case of magnetic field recordings, 
there are also advantages to having a static detector. Most importantly, if vector rather 
than scalar magnetometers are used (i.e. measuring individual directional components 
of the magnetic field), any variation in the orientation of the detector antenna 24 will 
15 again lead to significant spurious detector signals, since the magnetic detecting 
element will detect changes in the component of the geomagnetic field aligned with it. 
As a consequence of these effects, any translational movement of an electric field 
sensor comprising the detector antenna 24 or rotational movement of a magnetic field 
sensor comprising the detector antenna 24 will result in contamination of the detector 
20 signals by motionally induced noise. 

For these reasons it is desirable to cany out a controlled source 
electromagnetic survey to investigate or detect the presence of buried hydrocarbons 
using a combination of a mobile EM source operated just above the seafloor 6 and an 
array of electric and/or magnetic field sensors comprising the detector antenna 24 
25 placed statically on the seafloor 6 as indicatedin Figure 1 1 . At the ~end~of an EM 
survey, the instrument packages associated with the detectors are recovered using an 
acoustically actuated release mechanism to separate it from the ballast weight, 
allowing them to return to the sea surface for recovery, following standard 
oceanographic and marine geophysical practice. 
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Source and detector positions can be determined using standard long and/or 
short and/or ultra-short baseline acoustic navigation systems and processed to 
determine the separation between the source and detector. With traditional HED 
antennae, the relative orientation of a source and a detector must also be determined 
5 since HED detector signals depend critically on the relative geometry. This can be a 
substantial source of error when surveying with HED antennae. However, since the 
fields of both a VED antenna and VMD antenna are cylindrically symmetric around 
the source (see Figures 3a and 3b) only the source-detector separation is required 
meaning that geometric errors are reduced compared to those encountered using a 
1 0 traditional HED antenna. 

Furthermore, when using a towed HED antenna, both in-line and broadside 
geometries are required in order to unambiguously determine the subterranean strata 
resistivity structure. Although in-line data are straightforward to collect, ensuring the 
same coverage in the broadside data requires significantly more survey time since, as 
15 noted above, data can only be reliable collected with the HED antenna positioned at 
certain discrete positions (or survey locations) relative to a detector array and the HED 
antenna must also make multiple passes at each of these positions. An advantage of 
using a VED antenna and a VMD antenna is that detector signals resulting form the 
both the inductively coupled mode and the galvanically coupled mode can be 
20 collected continuously and simultaneously at each detector during towing, for 
example by driving the VED and VMD antennae at different frequencies. In some 
circumstances (for instance where the subterranean strata comprise substantially 
isotropic and homogeneous layers), it can be possible to separate VED and VMD 
detector signals collected simultaneously even with VED and VMD antennae being 
25 driven at the same frequency. For example, and considering only electric field detector 
signal components in a horizontal plane, this can be achieved by resolving the detector 
signal along radial and azimuthal directions. As noted above, the VED and the VMD 
antennae give rise to only radial and azimuthal electric fields respectively and so the 
azimuthal and radial components relate independently to the VED and VMD EM 
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signals. However, if the subterranean strata does not comprise substantially isotropic 
and homogeneous layers then refraction can lead to mixing of the signals from the 
VED and VMD antennae. 

It is noted that if the subterranean strata does comprise substantially isotropic 
5 and homogeneous layers, and in the case that the VED and the VMD antennae are 
driven at different frequencies, the resultant of the electric field amplitude measured in 
the two orthogonal horizontal directions of the detector antenna at each frequency 
reflects the radial and azimuthal components (and hence response of the subterranean 
strata to the VED and VMD EM signal respectively) regardless of the orientation of 

10 the detector antenna. However, in general the subterranean strata will not comprise 
substantially isotropic and homogeneous layers and knowledge of the detector antenna 
orientation, for instance measured by a compass, will allow more thorough data 
analysis to be performed. 

This means that a 2D section of the subterranean strata can be determined 

1 5 with relatively few detectors and EM source tows, so reducing the time taken, and 
hence cost, for an exploratory EM survey. If a full 3D map of the subterranean strata 
is required, the relative simplicity of the electric and magnetic field geometries 
associated with VED and VMD antennae means that an image of a target structure can 
be built up much more straightforwardly than using data from existing HED antenna 

20 systems 

It will be understood that whilst the above description describes a towed VED 
antenna, the method would also be applicable in a fixed installation. For example, the 
method could be used to monitor changes a hydrocarbon reservoir from which 
hydrocarbon is being drawn. In such cases it will be appropriate to employ one (or 
25 more)- VED- antenna(e) in fixed-positions relative to a detector array rather than 
undertake frequent towed surveys. The VED antenna(e) could be anchored to the 
seafloor or suspended from a oil-rig platform, for example. In other examples, the 
VED antenna(e) could be placed in a well or borehole, e.g. geotechnical borehole. The 
One (or more) VMD antenna(e) in fixed positions could also be used in conjunction 
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with the VED antenna(e). However, since the subterranean structures to which VMD 
EM signals are sensitive are unlikely to change with time, the results of an initial 
VMD EM survey (or MT survey) may be sufficient without needing to reproduce the 
VMD (or MT) part of the EM surveying. In the case of a producing oil field, the 
5 subterranean structures are likely to be well known already from prior geophysical 
surveys and drilling results. Indeed, prior geophysical and geological information from 
the oil field can be used to construct a background model without the need to collect 
VMD (or MT) data at any stage. 



10 The effects of the seawater surface 



A further consideration in designing an efficient EM survey (for instance when 
determining the range of source-detector separations R over which measurements can 
reliably be made) is the effects of the seawater surface on the propagation of the VED 
1 5 and/or VMD EM signals. 

In shallow water depths, it is possible for EM signals from an EM source to 
follow a propagation path upwards through the seawater to the surface; horizontally 
through the air; and back down through the seawater to the seafloor detector. This 
'airwave' component contains no information about the subterranean resistivity and so 
20 preferably should not dominate the detector signal. 

When using a HED antenna, the air wave component tends to dominate the 
detector signals in shallow water and at long source-detector separations. The higher 
the frequency of a drive signal applied to a HED antenna, the shorter the offset at 
which the airwave begins to dominate the detector signals. The effect of the airwave 
25 can be minimised by choosing appropriate transmission frequencies, and by targeting 
EM surveys on prospects in deep water and in which the target is at a relatively 
shallow depth below the seafloor [6]. 

The effect of the seawater surface is very different for a VED EM signal than 
for the signal from a HED antenna. A VED antenna induces predominantly vertical 
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current loops in the seawater. Because the vertical component of electric field is 
discontinuous at the sea surface (i.e. electric currents cannot flow across the interface), 
instead of propagating through the air, the EM energy carried in the EM signal is 
reflected back into the water. The effect of the seawater surface on a VED EM signal 
5 is therefore to introduce a surface-reflected signal component (as opposed to an 
airwave signal component) into the VED detector signal. Because the surface- 
reflected signal is rapidly attenuated by the conducting seawater, the effect of the 
seawater surface on a VED EM signal is more apparent at shorter source-detector 
separations than at longer source-detector separations. 

10 Figure 12a is a graph schematically showing model normalised vertical electric 

field component amplitudes, Ny 5 seen at a detector in response to a VED EM source 
signal as a function of separation, or range, R, between the VED antenna and the 
detector for a number of different overlying depths of seawater. The vertical electric 
field components are calculated for a uniform isotropic half-space subterranean strata 

15 configuration with resistivity 1 Qm, and seawater with resistivity 0.3 Qm, and a depth 
as indicated on each of the curves in the figure. In each case the electric field 
components are normalised by the modelled electric field components seen at a 
detector with the normalising model subterranean strata configuration described above 
(i.e. a uniform isotropic half-space subterranean strata configuration with resistivity 1 

20 Qm and seawater with resistivity 0.3 Qm and infinite depth). As before, the VED AC 
drive signal is a quasi-square wave at a frequency of 0.25 Hz and the electric fields ;are 
calculated per unit source electric dipole moment. The VED antenna is 100 m long 
and supported such that the lower electrode is 10 m above the seafloor. 

Figure 12b is a graph schematically showing model normalised radial electric 

25 field component amplitudes, N R , seen at a detector in response to a VED EM signal as 
a function of separation, or range, R, between the VED antenna and the detector for a 
number of different overlying depths of seawater. Figure 12b is otherwise similar to 
and will be understood from the description of Figure 12a above. 
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It is apparent from Figures 12a and 12b that the effects of a finite water depth 
on VED detector signals are most apparent in shallow water and at short source- 
detector separations. Larger values of Ny or Nr indicate that a larger fraction of the 
detector signal is due to the surface-reflected signal component. For example, at a 
5 detector-source separation of around 1000 m, and a water depth of 0.6 km, the radial 
component of the electric field is approximately double that seen at the same 
separation with an infinite depth of water (i.e. Nr = 2), the increase in detector signal 
being due to the surface-reflected signal component. 

As seen in Figures 5a, 5b and 7, the presence of a hydrocarbon layer in an 

10 otherwise homogeneous subterranean strata configuration becomes apparent at source- 
detector separations beyond around 4000 m and 6000 m for the radial and vertical 
electric field components seen at a detector in response to a VED EM source signal 
respectively. However, it can also be seen that the effects of finite seawater depth 
shown in Figures 12a and 12b are beginning to fall away beyond these detector-source 

15 separations (i.e. Ny and Nr are approaching unity), especially for seawater depths 
greater than around 0.6 km. This is in contrast to a HED antenna for which finite 
seawater depths can strongly affect detector signals at source-detector separations 
similar to those at which the detector signals are typically most sensitive to a 
hydrocarbon layer or reservoir. Accordingly, surveys which employ a VED antenna 

20 can be satisfactorily performed in shallower water than those employing conventional 
HED systems. 

For the VMD EM signal, the airwave component is present. It is similar to that 
observed with an HED, in that it dominates at long range and high frequency (for a 
given water depth) [15]. However since the response to a VMD is sensitive only to the 
25 background structure, and not to the hydrocarbon layer the VMD airwave component 
does not present a major problem. 



On the detection of the edges of a hydrocarbon layer by this technique 
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The model subterranean strata configurations used above to show how a 
combination of a VED and a VMD antenna can be used to detect a hydrocarbon layer 
and distinguish it from a general increase in resistivity with depth have all been ID 
models (i.e. horizontal layers of infinite horizontal extent). However, VED and VMD 
5 antennae EM surveying also works well in cases where a hydrocarbon layer comprises 
a reservoir of finite horizontal extent, and can also be used as a reliable means of 
locating the edges of an already identified reservoir structure. 

Figure 13a shows in schematic vertical section a fourth, or 3D model 
subterranean strata configuration. A VED antenna 22, a VMD antenna 23 and a 

10 detector antenna 24 are shown. The 3D model subterranean strata configuration 
includes a section of seafloor 6 beneath an infinite depth of seawater 4 with resistivity 
0.3 Qm. The strata beneath the seafloor 6 comprise a finite-extent hydrocarbon layer 
50 within an otherwise uniform background structure 52 of infinite horizontal and 
semi-infinite vertical extent. The uniform background structure has a resistivity of 2 

15 Qm. The finite-extent hydrocarbon layer has a vertical thickness of 100 m and an 8 x 
8 km square extent in a horizontal plane, its upper face is 1 km below the seafloor, it 
has a resistivity of 100 Qm and its centre is directly below the VED and VMD 
antennae. 

Figure 13b shows a schematic horizontal section through the centre of the 
20 finite-extent hydrocarbon layer 20 within the 3D model subterranean strata 
configuration shown in Figure 13a. The projected position of the VED and VMD 
antennae are marked by a cross 54. 

Figure 14 is a graph schematically showing the logarithm of the modelled 
radial and vertical electric field component amplitudes, Logio(E), seen by a detector in 
25 response to a VED EM signal as a function of separation, or range, R, between the 
VED antenna and the detector for various subterranean strata configurations. The 
separation between the source and the detector is made horizontally at the seafloor 
along a direction within the plane of Figure 13a and as denoted by the arrow marked R 
in Figure 13b. In this example, the VED AC drive signal is a quasi-square wave at a 
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frequency of 0.5 Hz and the electric fields are calculated per unit source electric 
dipole moment. The VED antenna is modelled as a point dipole source and supported 
10 m above the seafloor. The curves marked 3D-ver and 3D-rad respectively show the 
modelled radial and vertical electric field component amplitudes seen in response to 
5 the VED EM signal for the 3D model subterranean strata configuration shown in 
Figures 13a and 13b. For comparison, curves similar to those shown in Figure 5a (i.e. 
detailing the response of a normalising model subterranean strata configuration 
(dashed lines) and a ID hydrocarbon model subterranean strata configuration (solid 
lines)) are shown, but for consistency are computed for a 0.5 Hz VED AC drive signal 
10 and a background structure resistivity of 2 flm, as opposed to a 0.25 Hz VED AC 
drive signal and 1 Clm background structure resistivity. As in Figure 5a, these 
comparison curves are marked ver and rad for the vertical and radial electric field 
components, and shown as solid and dashed lines for the hydrocarbon and normalising 
model subterranean strata configurations respectively. 
15 The effects of the finite extent hydrocarbon layer on an otherwise uniform 

background are clear in Figure 14 from the difference between the curves marked 3D- 
ver and 3D-rad and the dashed-line curves marked ver and rad. The effects of the 
finite extent of the finite extent hydrocarbon layer are clear in Figure 14 from the 
difference between the curves marked 3D-ver and 3D-rad and the solid-line curves 
20 marked ver and rad. When both the source and detector are above the finite extent 
hydrocarbon layer (i.e. at source-detector separations R < 4 km) the curves 
representing the 3D model subterranean strata configuration (marked 3D-ver and 3D- 
rad) closely follow the curves representing the ID hydrocarbon model subterranean 
strata configuration (solid-lines marked ver and rad). When the detector is not above 
25 the finite extent hydrocarbon layer (i.e. at source-detector separations R > 4 km), the 
curves representing the 3D model subterranean strata configuration begin to approach 
those representing the normalising model subterranean strata configuration (dashed- 
lines marked ver and rad). 



WO 2004/008183 




PCT/GB2003/002379 



-43- 



Signal Requirements for Surveying 

The EM surveying methods described above can be readily applied to surveys 
for hydrocarbon resources on continental margins: 
5 The required EM source and detector characteristics in terms of power, signal 

to noise ratios and operating parameters will depend on application and can in general 
be met by existing instrumental technology, as can be seen by comparing the modelled 
electric field amplitude with recent publications from the academic sector. Seafloor 
electric field detectors can routinely measure electric field detector signals of 1(T U 
10 V/m [7, 14]. 

It can be seen from Figures 5a, 9a and 14 that at source-detector separations of 
around 8000 m, where the modelled VED detector signal is around 10* 15 V/m per unit 
source electric dipole moment, a detector signal of 10" n V/m can be obtained by 
driving the VED antenna with a VED AC drive signal capable of producing an 
15 electric dipole moment of around 10,000 Am. For the VED antenna described above, 
a suitable VED AC drive signal can be provided by a signal power supply unit and 
cycloconverter similar to those employed in a deep-towed active source instrument 
(DASI) currently in use [1, 8]. 

As can be seen from Figures 6a and 10a, the order of magnitude of the VMD 
20 detector signal electric field strengths are generally lower than the VED detector 
signal electric field strengths. Accordingly, a larger source dipole moment is required 
to provide detector signals which are detectable at appropriate source-detector 
separations. For example, at a source-detector separation of around 8000 m, the 
modelled VMD detector signal is around 10" 17 V/m per unit source magnetic dipole 
25 moment; Accordingly, to distinguish between different subFeiranean strata 
configurations, a VMD antenna magnetic dipole moment of 10 6 Am 2 is required. This 
could be accomplished, for example, with a 3m diameter loop VMD antenna 
comprising 1000 turns being driven by a VMD AC drive signal capable of providing a 
current of 150 A. The VMD antenna presents a substantially inductive load, however, 
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power supplies for driving such loads are readily available, for instance power 
supplies currently used for electric trams or train motors, or other heavy industrial 
machinery, which present a high inductance load could be used. 

5 Alternative Embodiment Based on VED Surveying Only 

According to an alternative embodiment, it is sufficient to collect VED data 
only. Instead of comparing the VED survey data set with a VMD (or MT) survey data 
set after normalisation of each set, it is possible to compare the VED data only with a 
10 suitable normalisation data set. In the above it has been assumed that the 
normalisation data sets do not have high integrity as regards resistivity profile 
especially, as will usually be the case when an area is poorly surveyed, for example 
when conducting initial exploration. However, in some cases an area will be very well 
characterised by previous surveying. For example, in a producing oilfield there is 
15 likely to be a wealth of existing seismic and well-log data. In this case, a background 
rock formation model can be created from the seismic data and then resistivities 
assigned to the various components in the rock structure using the resistivities 
obtained from well-log data. (If directly applicable well-log data is not available, it 
may be possible to estimate resistivity values by comparison with resistivity data from 
20 nearby wells in similar geological structures.) For these reasons, the present 
embodiment will be especially suited to applications in existing oilfields, such as 
monitoring long term depletion of reserves, as described above under the heading 
edge detection. 

Since the geological/geophysical data, such as seismic and well-log data, input 
into the rock formation model in this embodiment will be essentially blind to the 
presence of hydrocarbon, it may be thought of as performing an analogous role to the 
VMD or MT data in the previously embodiments. 

When monitoring depletion, it may be sufficient to directly compare VED 
survey data sets taken at different times, e.g. several weeks or months apart, without 
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use of normalisation. However, normalisation will generally be desirable, either using 
VMD, MT or rock formation model data as discussed above. 

Conclusions 

5 

The above description shows how controlled source electromagnetic surveying 
of subterranean strata using a VED antenna and a VMD antenna (or other source for 
inducing horizontal current flows) can allow the technique to be successfully applied 
10 to the problem of detecting hydrocarbon layers or reservoirs beneath the seafloor. In 
particular, we have identified that: 

A high quality survey can be carried out using a combination of a mobile EM 
source including a VED antenna and a VMD antenna operated just above the sea bed, 
and an array of detectors including electric and/or magnetic field sensing detector 
1 5 antenna placed statically on the seafloor. 

The effects of a finite seawater depth are less severe than seen when using a 
HED antenna. This allows a VED antenna to be operated in shallower seawater and in 
other circumstances where a HED antenna might not be useable. 

For purposes of data presentation and interpretation it is helpful to employ 
20 detector signal parameters that have been normalised with reference to an appropriate 
simplified model of the subterranean strata. 

In order to resolve the presence of a hydrocarbon layer, and to distinguish the 
effect of such a layer from other likely subterranean strata configurations, it is 
important to record detector signals seen in response to both inductively (i.e. 
25 horizontal current flow) and galvanically (i.e. vertical current flow) coupled modes of 
induced cuiTeht flowT 

3D modelling shows that with an appropriately positioned array of detectors, 
and with an appropriate source tow track, the method can yield valuable information 
about the limits of the areal extent of any subterranean hydrocarbon layer, as well as 
30 detecting its presence. 
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Using a VED antenna, the method can also be extended to EM sounding in 
wells, for use in 4D reservoir monitoring. This would also improve the coupling of the 
VED EM signals into the subterranean strata 

It will be understood that whilst the above description has concentrated on 
5 dipolar EM sources, similar methods employing other appropriately configured higher 
order multipoles, for example quadrupoles or octopoles could also be used. 

Finally it will also be understood that the invention is equally applicable to 
surveying of freshwater, for example large lakes, so that references to seafloor, 
seawater etc. should not be regarded as limiting. 
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CLAIMS 

1. An electromagnetic survey method for surveying an area that is thought or 
known to contain a subterranean hydrocarbon reservoir, comprising: 
5 obtaining a first survey data set from a vertical electric dipole (VED) source; 

and 

obtaining a second survey data set from a vertical magnetic dipole (VMD) 

source. 

10 2. An electromagnetic survey method according to claim 1, wherein the VED 
source and the VMD source are mounted on a common submersible vehicle. 

3. An electromagnetic survey method according to claim 1 or 2, wherein the first 
and second survey data sets are obtained together. 

15 

4. An electromagnetic survey method according to claim 1 or 2, wherein the first 
and second survey data sets are obtained separately. 

5. An electromagnetic survey method according to any one of the preceding 
20 claims, wherein the VED and VMD sources are operated at different frequencies. 

6. An electromagnetic survey method for surveying an area that is thought or 
known to contain a subterranean hydrocarbon reservoir, comprising: 

obtaining a first survey data set from a vertical electric dipole (VED) source; 

25 and 

obtaining a second survey data set from naturally occurring magnetotelluric 
(MT) electromagnetic fields. 
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7. A method of analysing results from an electromagnetic survey of an area that 
is thought or known to contain a subterranean hydrocarbon reservoir, comprising: 

providing a first survey data set obtained from a vertical electric dipole (VED) 

source; 

5 providing a second survey data set obtained from a vertical magnetic dipole 

(VMD) source; 

generating a first normalisation data set specific to the first survey data set; 
generating a second normalisation data set specific to the second survey data 

set* 

10 ' combining the first survey data set and first normalisation data set to obtain a 
first results data set that represents a difference between the first survey data set and 
the first normalisation data set; and 

combining the second survey data set and second normalisation data set to 
obtain a second results data set that represents a difference between the second survey 

1 5 data set and the second normalisation data set. 

8. A method of analysing results from an electromagnetic survey of an area that 
is thought or known to contain a subterranean hydrocarbon reservoir, comprising: 

providing a first survey data set obtained from a vertical electric dipole (VED) 

20 source; 

providing a second survey data set obtained from naturally occurring 
magnetotelluric (MT) electromagnetic fields; 

generating a first normalisation data set specific to the first survey data set; 
generating a second normalisation data set specific to the second survey data 

25 set; 

combining the first survey data set and first normalisation data set to obtain a 
first results data set that represents a difference between the first survey data set and 
the first normalisation data set; and 
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combining the second survey data set and second normalisation data set to 
obtain a second results data set that represents a difference between the second survey 
data set and the second normalisation data set. 

5 9. A method of analysing results from an electromagnetic survey according to 
claim 7 or 8, further comprising: 

normalising each of the first and second survey data sets relative to the 
respective first and second normalisation data sets or first and second functions 
specific to the first and second data sets respectively. 

10 

10. A method of analysing results from an electromagnetic survey according to 
any one of claims 7 to 9, wherein the first and second normalisation data sets or 
functions are calculated from a rock formation model. 



15 11. A method of analysing results from an electromagnetic survey according to 
any one of claims 7 to 9, wherein the first and second normalisation data sets or 
functions are calculated from the first and second survey data sets. 

12. A method of analysing results from an electromagnetic survey according to 
20 any one of claims 7 to 10, wherein the first results data set represents the difference 
between the first survey data set and the first normalisation data set as a function of 
position within the area, and the analysis of the first results data set includes 
identifying a location of a boundary of the subterranean hydrocarbon reservoir. 

25 13. A computer program product bearing machine readable instructions for 
implementing a method of analysing results from an electromagnetic survey according 
to any one of claims 7 to 12. 
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14. A computer apparatus loaded with machine readable instructions for 
implementing the method of analysing results from an electromagnetic survey 
according to any one of claims 7 to 12. 

5 15. A method of planning an electromagnetic survey of an area that is thought or 
known to contain a subterranean hydrocarbon reservoir, comprising: 

creating a model of the area to be surveyed including a seafloor, a rock 
formation containing a postulated hydrocarbon reservoir beneath the seafloor, and a 
body of water above the seafloor; 

10 setting values for water depth, depth below the seafloor of the postulated 

hydrocarbon reservoir, and resistivity structure of the rock formation; and 

performing a simulation of an electromagnetic survey in the model of the 
survey area by obtaining a first survey data set from a simulated vertical electric 
(VED) dipole source and a second survey data set from a simulated vertical magnetic 

1 5 dipole (VMD) source. 

16. A method of planning an electromagnetic survey of an area that is thought or 
known to contain a subterranean hydrocarbon reservoir, comprising: 

creating a model of the area to be surveyed including a seafloor, a rock 
20 formation containing a postulated hydrocarbon reservoir beneath the seafloor, and a 
body of water above the seafloor; 

setting values for water depth, depth below the seafloor of the postulated 
hydrocarbon reservoir, and resistivity structure of the rock formation; and 

performing a simulation of an electromagnetic survey in the model of the 
25 survey area by obtaining a first survey data set from a simulated vertical electric 
(VED) dipole source and a second survey data set from simulated magnetotelluric 
(MT) electromagnetic fields. 
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17. A method of planning an electromagnetic survey according to any one of 
claims 15 or 16, further comprising: 

repeating the simulation for a number of source-to-detector distances and 
frequencies in order to select optimum surveying conditions in terms of source-to- 
5 detector distance for probing the hydrocarbon reservoir. 

18. A method of planning an electromagnetic survey according to any one of 
claims 15 to 17, wherein the model includes a body of air above the body of water, 
and wherein the simulation takes account of signal propagation paths including the 

10 body of air when calculating the first and second survey data sets. 

19. A method of planning an electromagnetic survey according to any one of 
claims 15 to 17, further comprising: 

normalising each of the first and second survey data sets relative to respective 
15 first and second normalisation data sets or functions specific to the first and second 
survey data sets respectively. 



20. A computer program product bearing machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any one 

20 of claims 15 to 19. 

21. A computer apparatus loaded with machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any one 
of claims 15 to 19. 

25 

22. A submersible vehicle for subsea electromagnetic surveying, comprising: 
a vertical electric dipole (VED) source; and 

a vertical magnetic dipole (VMD) source. 
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23. A submersible vehicle according to claim 26, wherein the VED source and the 
VMD source comprise respective antennae mounted such that their dipole axes are 
aligned. 

5 24. A submersible vehicle according to claim 22 or 23, further comprising at least 
one waveform generator operable to transform a high voltage, low current AC drive 
signal received from an umbilical cable into a low voltage, high current AC drive 
signal to drive the VMD and VED antennae, 

10 25. A survey apparatus comprising: 

a submersible vehicle according to any one of claims 22 to 24; 

a signal power supply unit for generating a high voltage, low current signal 
power supply for the submersible vehicle; and 

an umbilical cable releasably connectable at ends thereof to the signal power 
1 5 supply unit and the submersible vehicle. 

26. A survey apparatus according to claim 25, further comprising a plurality of 
electromagnetic signal detectors. 

20 27. A surface vessel carrying a survey apparatus according to claim 25 or 26. 

28. A static platform carrying a survey apparatus according to claim 25 or 26. 

29. A well carrying a survey apparatus according to claim 25 or 26. 

25 

30. A borehole carrying a survey apparatus according to claim 25 or 26. 
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31. A method of analysing results from an electromagnetic survey of an area that 
is thought or known to contain a subterranean hydrocarbon reservoir, comprising: 

providing a survey data set obtained from a vertical electric dipole (VED) 

source; 

5 generating a normalisation data set specific to the survey data set, wherein the 

normalisation data set is calculated from a rock formation model; and 

combining the survey data set and normalisation data set to obtain a results 
data set that represents a difference between the survey data set and the normalisation 
data set. 

10 

32. The method of claim 31, wherein the rock formation model includes 
resistivity. 

33. The method of claim 31, wherein the rock formation model is derived from a 
1 5 combination of geological data and resistivity data. 

34. The method of claim 33, wherein the geological data is from seismological 
surveying. 

20 35. The method of claim 33 or 34, wherein the resistivity data is from well 
logging. 

36. The method of any one of claims 31 to 35, wherein the results data set 
represents the difference between the survey data set and the normalisation data set as 
25 a function of position within the area, and the analysis of the results data set includes 
identifying a location of a boundary of the subterranean hydrocarbon reservoir. 
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37. A computer program product bearing machine readable instructions for 
implementing a method of analysing results from an electromagnetic survey according 
to any one of claims 3 1 to 36. 



5 38. A computer apparatus loaded with machine readable instructions for 
implementing the method of analysing results from an electromagnetic survey 
according to any one of claims 3 1 to 36. 



39. A method of planning an electromagnetic survey of an area that is thought or 
10 known to contain a subterranean hydrocarbon reservoir, comprising: 

creating a model of the area to be surveyed including a seafloor, a rock 
formation containing a postulated hydrocarbon reservoir beneath the seafloor, and a 
body of water above the seafloor; 

setting values for water depth, depth below the seafloor of the postulated 
15 hydrocarbon reservoir, and resistivity structure of the rock formation; and 

performing a simulation of an electromagnetic survey in the model of the ' 
survey area by obtaining a survey data set from a simulated vertical electric (VED) 
dipole source; and 

adjusting the model to remove the postulated hydrocarbon reservoir and 
20 repeating the simulation to obtain a normalisation data set for comparison with the 
survey data set. 



40. A method of planning an electromagnetic survey according to claim 39, further 
comprising: 

repeating the simulation for a number of source-to-detector distances and 
frequencies in order to select optimum surveying conditions in terms of source-to- 
detector distance for probing the hydrocarbon reservoir. 
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41. A method of planning an electromagnetic survey according to claim 39 or 40, 
wherein the model includes a body of air above the body of water, and wherein the 
simulation takes account of signal propagation paths including the body of air when 
calculating the survey and normalisation data sets. 

42. A computer program product bearing machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any one 
of claims 39 to 41. 



10 43. A computer apparatus loaded with machine readable instructions for 
implementing the method of planning an electromagnetic survey according to any one 
of claims 39 to 41. 



44. A method of monitoring an area that contains a subterranean hydrocarbon 
1 5 reservoir by electromagnetic surveying, comprising: 

obtaining a first VED survey data set of the area; 

waiting a period of time; 

obtaining a second VED survey data set; and 

combining the first and second VED survey data sets to obtain a results data 
20 set that represents a difference between the first and second VED survey data sets, 
said difference being indicative of changes of the reservoir during said period of time. 

45. The method of claim 44, wherein the first and second VED survey data sets 
are also combined with a normalisation data set. 

25 

46. The method of claim 45, wherein the normalisation data set is obtained from a 
survey of the area by a vertical magnetic dipole (VMD) source 
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47. The method of claim 45, wherein the normalisation data set is obtained from a 
survey of the area by naturally occurring magnetotelluric (MT) electromagnetic fields. 



48. The method of claim 45, wherein the normalisation data set is calculated from 
5 a rock formation model. 

49. The method of any one of claims 44 to 48, wherein the period of time is at 
least a day, a week, or a month. 

10 50. The method of any one of claims 44 to 49, further comprising: analysing the 
results data set to identify changes in boundaries of the reservoir during the period of 
time. 

51. A computer program product bearing machine readable instructions for 
1 5 implementing the method of any one of claims 44 to 50. 

52. A computer apparatus loaded with machine readable instructions for 
implementing the method of any one of claims 44 to 50. 



20 



